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Abstract 

Modeling the intermittent behavior of turbulent energy dissipation processes both in space and 
time is often a relevant problem when dealing with phenomena occurring in high Reynolds number 
flows, especially in astrophysical and space fluids. In this paper, a dynamical model is proposed to 
describe the spatio-temporal intermittency of energy dissipation rate in a turbulent system. This is 
done by using a shell model to simulate the turbulent cascade and introducing some heuristic rules, 
partly inspired by the well known p-model, to construct a spatial structure of the energy dissipation 
rate. In order to validate the model and to study its spatially intermittency properties, a series 
, of numerical simulations have been performed. These show that the level of spatial intermittency 

O ■ 

of the system can be simply tuned by varying a single parameter of the model and that scaling 
laws in agreement with those obtained from experiments on fully turbulent hydrodynamic flows 
£> . can be recovered. It is finally suggested that the model could represent a useful tool to simulate 



> 

in 



X 



the spatio-temporal intermittency of turbulent energy dissipation in those high Reynolds number 
astrophysical fluids where impulsive energy release processes can be associated to the dynamics of 
the turbulent cascade. 

PACS numbers: 47.27.Eq, 02.50.Ey, 47.53. +n 
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I. INTRODUCTION 



The dynamics of fluids and plasmas, both in laboratory experiments and in astrophysical 
or geophysical systems, is very often characterized by the presence of turbulent motions |l|, Q] ■ 
In several contexts of astrophysics and space physics, it is extremely important to model 
some of the effects related to the turbulent dynamics. In particular, describing in a proper 
way the spatio-temporal intermittency of the turbulent energy dissipation process is one of 
the basic ingredients for the study of several astrophysical systems. As relevant examples, 



we can consider the active regions of the solar corona |3j, |4j, the interstellar medium [5], and 
accretion disks 



6]. 

Intermittency is one of the most investigated problems in the field of fully developed 
turbulence (see jl| and references therein). Among the many approaches used for the study 
of intermittency in turbulence, here we want to briefly recall only some of them, which are 
related to the work presented in this paper. A number of random cascade models (see e.g. [?], 

[2I, Eo]) were initially proposed to reproduce the observed intermittency corrections (see 
e.g. pjj) to the scaling laws of the classical Kolmogorov theory of turbulence jl2j . Another 
interesting approach to the modeling of intermittency of the turbulent energy cascade is 
based on the use of dynamical deterministic models known as shell models (see the reviews 



by Bohr et al. 



131 ]. Giuliani [lj], Biferale 151]) . More recently, several relevant developments 



have led to a beginning of a deeper understanding of the intermittency phenomenon. To 



mention but a few: the role of Lagrangian conservation laws 16j and nonlocal interactions 



171 ] on intermittency, and the introduction of new multifractal approaches for the description 



of velocity increments statistics 



18|- 



Besides these important theoretical advances, there are several more specific situations 
where a simple dynamical system modeling of the intermittency in the turbulent cascade 
can be extremely helpful. This can be the case of astrophysical and space fluids, where, due 
to the extremely large Reynolds numbers, dynamical models which are able to simulate the 
turbulent cascade and the related energy dissipation processes in Reynolds number regimes 
which are not far from the real ones (at least with respect to direct numerical simulations) can 
represent an essential ingredient for the modeling of such physical systems. An example is 
given by the recent applications of shell models to the description of the statistical properties 



of solar flares 



19 



201 ] . and to the nanoflares occurring in solar coronal loops 2jJ]. In this 
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framework, it is worth to point out that shell models provide only a temporal description 
of the intermittency properties since they lack any spatial information. The possibility to 
have a dynamical "shell-like" model capable of reproducing some intermittency properties 
both in space and time would thus be attractive. 

For the reasons explained above, in this work we propose a simple method to model the 
intermittent character of energy dissipation in a turbulent system in both space and time, 
by using a shell model together with some rules inspired to some extent by the well known 
turbulence p-model lp| . The paper is organized as follows. In Section [Til we recall the main 
ideas concerning shell models and the p-model, in Section II III we give a description of the 
proposed method, in Section HVl we provide some details about the numerical procedure, in 
Section |V] we show the results of the spatially intermittency analysis performed in order to 
validate the proposed model, while the conclusions are drawn in Section I VI I 



II. BACKGROUND 



Shell model 



Shell models were introduced in 70's by Obukhov 



221 ]. Gledzer [23J, and Desnyansky and 



Novikov [24J in the context of hydrodynamic turbulence and, since then, used extensively 



261 ] ) and magnetohydrodynamics (MHD) (see e.g. 



both in hydrodynamics (see e.g. [25|, 

2?1 [28, 29, 30]). They are based on a set of coupled nonlinear ordinary differential equations 
which describe the dynamics of the turbulent energy cascade in the wave-vector space. The 



dynamical and statistical behavior of shell models have been investigated in detail in many 
works (see [l3, Q, [ijj] and references therein) and it has been shown that they are able to 
describe several properties of the turbulent energy cascade process. The main advantage of 
shell models is that they can be investigated through numerical simulations at high Reynolds 
numbers much more easily than Navier-Stokes (N-S) or MHD equations, due to the reduced 
number of degrees of freedom. On the other hand, an obvious minus of these models is the 
absence of any information about spatial structures. 

Shell models are built up by dividing the wavevector space in a discrete number of shells 
of radius k n = ko\ n , with A > 1 fixing the shell logarithmic spacing (usually A = 2) and 
n = 1,2, ...,N. Each shell is associated with a dynamical complex variable u n (t) which 
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represents the time evolution of velocity fluctuations at scale i n ~ k^ 1 . The evolution 
equations for the variables u n (t) are written by introducing nonlinear terms in the form of 
quadratic couplings between neighbouring shells. The coupling coefficients are chosen to 
satisfy scale invariance and the conservation of the ideal invariants, as for example the total 
energy and the kinetic helicity for N-S equations. 

For this work we use the hydrodynamic shell model proposed by L'vov et al. jsi], also 
known as Sabra model. The evolution equations of the shell variables are 

du n _ ( * i * , i \ 

—j— — %K n yU n+ 2U n+1 — jU n+ \U n _ l + g1t n -lW n -2 J 

-vk 2 n u n + /„ , (1) 

where v represents the kinematic viscosity, and /„ is a forcing term usually acting on some 
low wavenumber shells. From Eq. (pQ), we can derive the evolution equation for the n-th 
shell kinetic energy E n = (u n ii*)/2: 

dt 

where 



TI n - 2uk 2 n E n + 3fJ«/ n ) , (2) 



n n = k n Q(u* n u n+2 u* n+1 + \u n u* n+l u n -i 

+ \u* n U n - 1 U n -2) (3) 

gives the kinetic energy flux through the n-th shell (the symbols 9ft and S denote the real 
and imaginary parts respectively of a complex number). 
The total energy dissipation rate e(t) can be defined as 

N 

e(t)=uJ2kl\u n \ 2 . (4) 



n=l 



B. The p-model 

The p-model has been designed to describe the observed multifractal behavior of the 
energy dissipation rate in fully turbulent flows [10]. Without loss of generality, in this 
paper we consider, for simplicity, a one-dimensional spatial domain. In this case, the total 
dissipation rate e r (x) in the segment [x,x + r] is equal to £ r (x) = J^ +r e(x)dx, e(x) being 
the energy dissipation rate in the x position. In the p-model, an interval of size r breaks 



down into two sub intervals of size r/2, and the energy flux to these smaller eddies proceeds 
unequally. A fraction p (with 0.5 <p < 1) of the dissipation contained in the parent interval 
is distributed equally on one of the two subintervals (left or right with equal probability), 
and the remaining (1 — p) fraction on the other sub interval. This process starts from the 
integral scale L (where we have only one interval) and is repeted until segments of size r] 
(corresponding to the dissipative scale) are created. It has been shown in Ref. [h]] that 
using p = 0.7 the multifractal spectrum of the synthetic dissipation signal obtained through 
the p-model reproduces extremely well the results of experiments. 



III. THE MODEL 



In a few words, the basic idea of the method proposed here consists in using a shell 
model to describe the dynamics of the turbulent cascade process and in providing a spatial 
structure to the energy dissipation using some rules, which partly recall the p-model, to 
distribute in space the energy fluxes given by the shell model. 

We consider a one-dimensional spatial domain whose size is denoted by L. As a base for 
the construction of the spatial energy structure we use a hierarchy of N scales £ n = 2 1 ~ n L 
(n = 1,2, ...,N). For each scale n we can define a set of 2™ _1 disjoint segments of size £ n 
which cover the spatial domain. Let us note that this is the same hierarchical structure 
as in the p-model, that is, each segment at the scale n can be considered as parent of two 
corresponding segments at the scale n + 1 which have half the size and cover the same 
subinterval as the parent. 

Let us now suppose that the N scales of this hierarchy are associated with the N shells 
of the shell model, and that l n = l/k n . At each time step £j (i=l,2,...) of the numerical 
solution of the shell model equations, we can calculate the energy increment AE n (ti) of the 
n-th shell as 

AE n (ti) = E n (ti) - £ n (ti_ a ) . (5) 

These increments are used to construct a spatial energy structure which evolves in time 
parallely to the shell model as explained below. 

The increments AE n (ti) are distributed over the spatial grid of the corresponding scale 
using the following criteria: 
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• If AE n (ti) > we first divide the energy AE n (U) among the segments at the scale n — 1 
(which thus play the role of parent segments) proportionally to the energy contained in 
them at the scale n — 1. The energy increment thus obtained for each parent segment 
is then transferred to the corresponding two daughter segments in the same way as in 
the p-model, that is, adding a fraction p of the increment (0.5 < p < 1) to one of the 
daughters (left or right with equal probability) and the remaining (1 — p) fraction to 
the other daughter. 

• If AE n (ti) < the increment is subtracted from the energy at the scale n in such a 
way that each segment is depleted by a fraction of AE n (ti) proportional to the energy 
content of the segment itself at the same scale. 

As a result of the procedure described above, to each one of the 2 n_1 segments which 
cover the domain at the scale n is attributed a kinetic energy E\ n \t) (where / = 1, ...,2™ _1 
is the index denoting the segments at the scale n). The total energy at the scale n equals 
the kinetic energy of the corresponding shell, that is, ^^(0 = \ u n{t)\ 2 

In order to have an evolution of the spatial energy structure, the spatial distribution 
of the p and (1 — p) values is changed during the time evolution. Two different methods, 
described below, were used to perform this changes in time and we denote by Model A and 
Model B the two versions of the model corresponding respectively to these two methods. 

Model A. In the first version of the model the changes of the probabilities in space are 
done at the same time instant for all the segments at a given scale n. Let us suppose that 
the last change in the p distribution at the scale n occurs at the time step . At each time 
step ti we calculate the instantaneous eddy turnover time T^ a ~ l \ti) = [k n ^iU n -i(ti)Y l for 
the scale n — 1 and compare it to the time elapsed from the last change At^ — U — t^K 
If At^ > T^™ -1 ), the p spatial distribution at the scale n is redrawn. This procedure is 
followed for each n > 1. 

Model B. In the second version the spatial changes of the probabilities are done inde- 
pendently for the different segments at a given scale n. Indeed one of such changes always 
involves a couple of segments / and / + 1 (where / is an odd integer number), because if the 
probability is p for the segment / it must be 1— p for the segment l+l and viceversa. As a con- 
sequence, if we denote by t^ the last change for the segment I at the scale n, we have that 
^(n,o _ ^(n,i+i) a COU pj e Q f neighbouring segments with / odd integer. At each time step U, 
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we compute the local instantaneous eddy turnover time for the corresponding father segment 
at the scale n - 1 (whose index is l f = i±l), that is, r^\u) = [k n ^2E^\t^ 
and compare it to the time elapsed from the last change At^'Q = At*™'^ 1 ) = tj — t^ ,l \ If 
' s iti) the probability values are redrawn for the segments / and I + 1. This 
procedure is followed for each n > 1. 

The two procedures described above aim to describe phenomenologically the correlations 
arising in the cascade due to the scaling of the eddy turnover times. The main difference 
between them consists in the fact that in Model B we try to take into account also the local 
dynamics of the turbulent cascade on the spatial domain. 

In both the versions of the model, summing the contribution coming from all the scales, 
we obtain the spatial shape of the energy density w(x,t) as 

™0M) = E-7^' (6) 

n=l ln 

where x denotes the grid position, corresponding to the smallest scale grid spacing, and 
ln(x) = Int[(x — I) /i n ] + 1 (Int denotes the integer part of a real number). We can now 
define also an energy dissipation rate which depends on the spatial coordinate as 

e(x,t) = V ^k>-±@-. (7) 

n=l l n 

IV. NUMERICAL PROCEDURE 

The shell model equations Eqs. ([T]) have been numerically solved using a 4-th order 
Runge-Kutta scheme. The parameters used in the simulations are N — 15, v — 10~ 4 , 
ko = 2 -5 , and A = 2. We used an external forcing term applied on the 3-rd and 4-th shell 
given by fs — f^ — 0.1(1 With these parameters the Reynolds number is Re ~ 10 5 and 
the large scale eddy turnover time r e ~ 40. A sample of the total energy dissipation rate 
e(t) given by the shell model is shown in Fig. [TJ 

The number of grid points in the 1-D spatial domain is 2 14 . Another free parameter of 
the model is the p value used to construct the spatial structure of the energy dissipation. It 
should be clear to the reader that it is not necessary to use the value p = 0.7 which allows 
the p-model to reproduce the multifractal structure of the energy dissipation rate observed 
in experiments of fully developed turbulence }lo| . This is why the models proposed here, 
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FIG. 1: A sample of the total energy dissipation rate given by the shell model. 



although inspired to some extent by the p-model, are substantially different from it, being 
characterized also by the dynamics provided by the shell model and by the evolution of the 
p spatial distribution. The multifractal properties of the spatial energy dissipation given by 
the model change in time. More precise indications about the values to be attributed to 
the p parameter can be found from the application of the model to well defined physical 
situations, a question which we plan to further investigate in the future. For this work, we 
have performed numerical simulations using the values of p — 0.7,0.8,0.9,1 for Model A, 
and p = 0.7, 0.8 for Model B. We would like to remark that changing the spatial distribution 
of the p and 1 — p values according to the instantaneous eddy turnover times at each scale, 
as described in the previous section, is necessary for describing the intermittent behavior of 
the energy dissipation. As a confirmation, we performed a simple test in which we modified 
the spatial distribution of the probabilities at each time step and obtained that in this case 
the energy dissipation rate is not intermittent for any value of p. 



V. ANALYSIS OF THE SPATIAL INTERMITTENCY 

In order to validate the proposed models, we performed an analysis of the spatial inter- 
mittency properties of the energy dissipation rate for different values of the p parameter. 
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This has been done both by calculatin g th e kurtosis of e(x) and by looking for the presence 
of multifractal scaling laws in e(x) 32|, l33l. l34| . In this way, some comparisons to the scaling 



exponents of structure functions found in turbulence experiments can also be made. 
The starting point for multifractal analysis is the definition of the probability measure 

*(0 = *g , (8) 

where 



Xi(r) = / e(x)dx . (9) 

JSi(r) 

S{(r) represents a hierarchy of disjoint segments of size r covering the domain L. We can 
calculate the so called generalized dimensions D q |35| by looking at the scalings of the g-th 
order moments of fii(r) vs. r: 

(A* 8 > = E^(r)~r^«. (10) 

i 

The largest values of q amplify the contribution given to (fi q ) by the most intermittent 
regions of the measure, while for small values of q the major contribution is due to the most 
regular regions. If the probability measure is globally self-similar (i.e. non intermittent), D q 
is constant and it corresponds to the fractal dimension of the measure. Conversely, if D q is 
not constant, the scaling laws are said to be anomalous and the measure can be described as 



a multi 



of q [33] 



ractal object. In this case, it can also be shown that D q is a nonincreasing function 



The generalized dimensions D q can also be related to the scaling exponents ( q of the 
velocity structure functions, which are measured in fluid flows and represent the benchmark 
for the nonlinear energy cascade modeling. These exponents are defined by 

(Svf) = ([v(x + r) -v(x)] 9 } ~r c * . (11) 



It can be shown 



361 ] that 



C 3 = I + (I - l) (D q/ s - d) , (12) 



where d represents the topological dimension of the support, in our case d — 1. 

We show in the next two subsection the results obtained for Model A and Model B 
respectively. 
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A. Model A 



In Fig [2] the space-time structure of the energy dissipation rate e(x,t), calculated accord- 
ing to Eq. (j7|), is shown for p = 0.8 and p — 1. For the sake of clarity the grey levels refer 
to the logarithm of e(x,t). It can be seen that e(x,t) becomes more and more fragmented 
in space as p increases, as one would expect. 

To give a first indication on the intermittency properties of the spatial energy dissipation 
e(x), the time evolution of the kurtosis of e{x) for the four different values of p used is 
shown in Fig. [31 The increase of the typical values of the kurtosis confirms that the level of 
intermittency is significantly enhanced as p increases from 0.7 up to 1. In Fig. [4] the spatial 
structure of the energy dissipation rate e(x) is shown for four fixed time instants (one for 
each of the different p values used) where the kurtosis shown in Fig. [3] displays a peak. The 
four time instants chosen are t = 73 for p = 0.7, t = 133 for p = 0.8, t = 80 for p = 0.9, 
and t = 51 for p = 1. This figure shows that at some positions the quantity e(x) shows very 
strong bursts, which appear to become stronger and stronger as p increases. 

The multifractal analysis was performed by calculating the moments (fi q ) given in Eq. 
(fIU|) for —5 < q < 5 at 100 different time instants, namely t = 41, 42, 140. A good scaling 
region extending almost over the whole r range was found for all the g's. The generalized 
dimensions D q were calculated as averages over all the time instants considered. The plot 
of D q vs. q obtained for the four values of p used is shown in Fig. [51 It can be seen that the 
spectrum of D q values becomes wider as p increases. This is a consequence of the enhanced 
level of intermittency which produces stronger and more localized dissipation bursts for 
larger values of p. 

The scaling exponents ( q as estimated from the energy dissipation scalings in Model A for 
the four values of p considered are reported in Table ID As a comparison, the ( q exponents 
computed from a wind tunnel experiment [37J] are also shown. It can be seen that Model A 
gives scaling exponents in agreement with experiments (within the experimental error) for 
p > 0.9. 
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FIG. 2: Space-time structure of the energy dissipation rate e(x,t) in Model A for p = 0.8 (top 
panel) and p = 1 (bottom panel). The grey levels refer to the logarithm of e(x,t). 
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FIG. 3: Time evolution of the kurtosis of the spatial energy dissipation rate e(x) in Model A for 
the four different values of p used: from top to bottom p = 0.7, p = 0.8, p = 0.9, and p = 1 
respectively. 
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FIG. 4: Energy dissipation rate on the spatial domain in Model A: from top to bottom, 
and t = 73,p = 0.8 and t = 133, p = 0.9 and t = 80, p = 1 and t = 51. 
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FIG. 5: Generalized dimensions D q of the spatial energy dissipation rates in Model A for p = 0.7 
(Dotted line), p = 0.8 (Dash- Dotted), p = 0.9 (Dashed line), and p = 1.0 (Solid line). The D q 
values were calculated as time averages over all the time instants considered (see text). Error bars, 
representing standard deviation errors, are also reported for some values of q. 

B. Model B 

The space-time structure of the energy dissipation rate e(x,t) for p = 0.7 and p = 0.8 is 
shown in Fig. [6j A larger spatial fragmentation of e(x, t) for the larger p is clearly observed. 

The time evolution of the kurtosis of e(x) for p = 0.7 and p = 0.8 is shown in Fig. fTJ It 
can be seen that the typical values of the kurtosis increase going from p = 0.7 top = 0.8 
as it could be expected. Comparing Fig. [TJto Fig. Owe can notice that, for the same p, 
Model B is much more intermittent than Model A. For instance, using p = 0.8 Model B 
roughly reaches the same level of intermittency as Model A with p — 1. In Fig. [8] we show 
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FIG. 6: Space-time structure of the energy dissipation rate e(x, t) in Model B for p = 0.7 (top 
panel) and for p = 0.8 (bottom panel). The grey levels refer to the logarithm of e(x,t). 
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TABLE I: Scaling exponents Q q for the velocity structure functions as estimated from the energy 
dissipation scalings in Model A for the four values of p considered. In the last column we report 
the velocity structure functions exponents computed from a wind tunnel experiment [37]. 



q 


p=0.7 


p=0.8 


p=0.9 


p=1.0 


wind tunnel 


1 


0.337 ±0.001 


0.342 ± 0.002 


0.350 ± 0.005 


0.362 ± 0.008 


0.37 ±0.01 


2 


0.670 ± 0.001 


0.675 ± 0.002 


0.683 ± 0.005 


0.694 ± 0.007 


0.70 ±0.01 


3 


1.00 


1.00 


1.00 


1.00 


1.00 


4 


1.326 ±0.002 


1.317 ±0.004 


1.304 ±0.008 


1.285 ±0.012 


1.28 ±0.02 


5 


1.650 ±0.005 


1.628 ±0.011 


1.596 ±0.020 


1.554 ±0.027 


1.54 ±0.03 



the spatial structure of the energy dissipation rate e(x) at two time instants where a peak 
in the kurtosis is found, that is, t = 58 for p = 0.7 and t = 66 for p = 0.8. From this figure 
it is clear that also for Model B e(x) shows very strong intermittency bursts in space. 

Also for Model B we investigated the spatial intermittency properties of the energy dissi- 
pation rate through the multifractal analysis described previously. The moments {n q ) given 
in Eq. ( JTU1) were calculated also in this case for —5 < q < 5 at the 100 time instants 
t = 41, 42, 140. Good scalings were found for all the q values and the dimensions D q were 
obtained as averages over all the time instants considered. Fig. [9] shows the plot of D q vs. 
q for p = 0.7 and p = 0.8. As expected, the D q spectrum is wider for p = 0.8. Moreover we 
can observe that the D q curve for p = 0.8 is very close to the one obtained in Model A for 
p — 1. 

The scaling exponents ( q of the velocity structure functions in Model B as estimated from 
Eq. (j!2p are reported in Table [Til A good agreement with the scaling exponents found in 
the wind tunnel experiment analyzed in Ref. 37J is obtained for p = 0.8. 

The fact that for a given p a larger intermittency is found in Model B than in Model 
A is clearly a consequence of the fact that in the procedure used in Model B for the time 
evolution of the spatial distribution of p and 1 — p we consider also the local dynamics of 
the energy cascade as pointed out in Section [TLT1 
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FIG. 7: Time evolution of the kurtosis of the spatial energy dissipation rate e(x) in Model B for 
p = 0.7 (top panel) and p = 0.8 (bottom panel). 

VI. CONCLUSIONS 



A problem which often arise when studying turbulent phenomena occurring in astrophys- 
ical and space fluids is the description of the intermittency of the energy dissipation process 
from a spatio-temporal point of view. Due to the huge Reynolds numbers occurring typically 
in these situations a dynamical system modeling of space-time intermittency can represent 
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FIG. 8: Energy dissipation rate on the spatial domain in Model B for p = 0.7 and t = 58 (top 
panel), p = 0.8 and t = 66 (bottom panel). 



an important ingredient for the characterization of such systems. 

In this paper we propose a method to model the main intermittency features of energy 
dissipation in a turbulent system both in space and time. This is done by using a turbulence 
shell model and introducing some heuristic rules, partly inspired by the well known cascade 
p-model, to construct a spatial structure for the energy dissipation rate. 

To the aim of validating the model, we performed a series of numerical simulations to 
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FIG. 9: Generalized dimensions D q of the spatial energy dissipation rates in Model B for p = 0.7 
(Dashed line) and p = 0.8 (Solid line). The D q values were calculated as time averages over all the 
time instants considered (see text). Error bars, representing standard deviation errors, are also 
reported for some values of q. The D q curve obtained in Model A for p = 1 is also shown as a 
comparison (dotted line). 

study the spatial intermittency properties of the energy dissipation rate for different values 
of the free parameter p. The results show that the level of spatial intermittency of the 
system can be simply tuned, in both the proposed versions of the model, by changing 
the value of p. The spatial intermittency of the system is enhanced by increasing the p 
parameter. Scaling laws in agreement with those obtained in experiments involving fully 
turbulent hydrodynamic flows are recovered in Model A for 0.9 < p < 1, and in Model B 
for p ~ 0.8. 

The results of this work open the way to applications of the proposed model to different 
physical situations. In our opinion this model could represent a useful tool to simulate 
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TABLE II: Scaling exponents ( q for the velocity structure functions as estimated from the energy 
dissipation scalings in Model B for p = 0.7 and p = 0.8. In the last column we report the velocity 
structure functions exponents computed from a wind tunnel experiment \st\ - 



Q 


p=0.7 


p=0.8 


wind tunnel 


1 


0.348 ± 0.004 


0.369 ± 0.009 


0.37 ±0.01 


2 


0.681 ±0.004 


0.700 ± 0.008 


0.70 ±0.01 


3 


1.00 


1.00 


1.00 


4 


1.307 ± 0.006 


1.277 ±0.014 


1.28 ±0.02 


5 


1.602 ±0.015 


1.536 ±0.031 


1.54 ±0.03 



the spatio-temporal intermittency of turbulent energy dissipation in those high Reynolds 
number astrophysical fluids where impulsive energy release processes can be associated to 
the dynamics of the turbulent cascade. To give just some examples, such a modeling could 
be interesting for studying the role of intermittent energy dissipation in the active regions of 
the solar corona, in the interstellar medium, and in accretion disks. We plan to investigate 
some of these problems in future studies. 
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